Pseudogap Kondo Physics from Charge Fluctuations in a Quantum Dot 
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We consider charge fluctuations in a quantum dot coupled to an interacting one-dimensional 
electron liquid. We find the behavior of this system to be similar to the multichannel pseudogap 
Kondo model. By tuning the coupling between the dot and the one-dimensional electron liquid, 
one can access the quantum critical point and the various fixed points which arise. The differential 
capacitance is computed and is shown to contain detailed information about the system. 



Nanotechnology has been the source of a renewed in- 
terest in the Kondo effect. The incredible progress in 
miniaturizing solid state devices has made it possible 
to fabricate small metallic islands (i.e. quantum dots) 
by confining electrons in a two-dimensional electron gas. 
Quantum dots provide a highly controllable environment 
to study Kondo physics, and allow for many aspects of 
the Kondo effect to be probed. In this work, we sug- 
gest that a quantum dot coupled to an interacting one- 
dimensional electron liquid (i.e. a Luttinger liquid) could 
provide a controlled environment to observe pseudogap 
Kondo physics. 

The pseudogap Kondo model was first considered in 
Ref. ^. In this model, a magnetic impurity is coupled 
to a sea of conduction electrons with a density of states 
vanishing at the Fermi energy with power-law behavior 



Hint = J T ■ 



-</4(0)V s (0) with p(E)=p \E\ 



One of the most interesting features of this model is that, 
for antiferromagnetic coupling (J > 0), there is an un- 
stable intermediate coupling fixed point occurring when 
J = J c . For J < J c , the impurity spin is unscreened at 
low energies; for J > J c the impurity spin is screened. 
At J = J c , the impurity spin exhibits quantum critical 
fluctuations. It is worth mentioning that this model has 
attracted attention recently due to its potential relevance 
to various correlated electron systems. In particular, this 
model has been argued to describe impurities in high-T c 
cuprate superconductors, fl Moreover, the critical be- 




I 

l 



FIG. 1: Setup: A quantum dot coupled to an interacting one- 
dimensional electron liquid. The number of electrons on the 
dot is controlled by the gate voltage V g . A voltage applied to 
the auxiliary gates, V a , controls the coupling between the dot 
and the one-dimensional electron liquid. 



havior occuring when J = J c may be relevant to the 
behavior seen in heavy fermion materials.^, 5] 

The setup we consider is shown in Fig. 1. A large 
quantum dot is coupled to a reservoir, consisting of an 
interacting one-dimensional electron liquid. The dot is 
capacitively coupled to a gate; the gate voltage V g con- 
trols the number of electrons on the dot. The coupling 
between the dot and the reservoir is controlled by a volt- 
age V a applied to the auxiliary gates. To model the dot, 
we assume the level spacing of the dot is much smaller 
than any energy scale in the problem; we approximate 
the spectrum of the dot by a single particle continuum. 
Moreover, we consider the case where the reservoir is cou- 
pled to the dot via a point contact. The Hamiltonian for 
the dot has the form f/dot = flo(dot) + flint- flo(dot) 
describes the single particle energy levels of the dot; H mt 
describes the charging energy of the dot, as well as the 
coupling to the reservoir 



(i) 



In Eq. [l], ip2,s (ipi,s) destroys an electron in the dot (reser- 
voir); N is the number operator of the dot; N is the 
average number of electrons on the dot, which is propor- 
tional to V g \ E c is the charging energy of the dot; t is 
the tunneling matrix element between the dot and the 
reservoir, which is controlled by V a . For generic values of 
N, it costs a finite energy to put an extra electron on the 
dot; for temperatures sufficiently less that E c , Coulomb 
blockade develops and the number of electrons on the dot 
becomes quantized. However, for N = n + 1/2 the ener- 
gies of the n-electron and (n-fl)-electron states are equal, 
and the charging energy vanishes. Therefore, quantum 
fluctuations between the dot and the reservoir become 
important. 

In the following, we assume t is small and we focus on 
the regime N m n + 1/2. For energies sufficiently less 
than Ec, the physics will be dominated by the states n 
and (n + 1). Hence, we can project out all other states 
and restrict ourselves to this subspace. By considering 
these states as the two states of a pseudospin t z = ±1/2 
and writing N = (n + 1/2) + t 2 , Hi n t takes the form[|| 



flh 
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m js (0) + h.c. 



hr z 



(2) 



where h = E c [N - (n + 1/2)]. Eq. | is a Kondo Hamil- 
tonian with anisotropic couplings. Whereas the Kondo 
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effect usually involves a magnetic impurity, it arises in 
this system due to charge fluctuations. 

Recently, it was argued that, besides the Kondo 
physics of Eq. 0, other types of behavior are possible. 
In particular, by performing a variational calculation, 
these authors identified that quantum fluctuations might 
give rise to tricritical Ising behavior. However, it remains 
to be seen whether these results will be confirmed nu- 
merically or experimentally. In this work, we focus on 
regimes where the system is far from the potential tri- 
critical point, so that the Kondo physics dominates. 

Being interested in the low energy properties of the 
system, we expand the electron operator in the reservoir 
in terms of right and left movers 



ikpx 



1pR,l,s( X ) + ' 



-ikpx 



where kp is the Fermi wavevector, and V>ii,i.s and ipL.i.s 
are the (slowly varying) right and left moving fermion op- 
erators. Moreover, upon expanding the electron operator 
in the dot in harmonics centered about the point con- 
tact, the reservoir couples to only a single harmonic. |8| 
Focussing on that single harmonic, we can write an ef- 
fective one-dimensional model for the dot.J^] In what 
follows, we will make extensive use of the boson rep- 
resentation. To do so, the electron operator is writ- 
ten as ipR/L.i.s ~ e ±lV ^^ R / t - < >* where the chiral fields, 
(f>R,i,a and 4>L,i,s, are related to the usual Bose field 
ef>i <a and its dual field 9 ilS by (j) ilS = (f>R,i )S + ef>L,i, a and 
Oi.s = <f>R,i,s ~~ 4>L,i,s- It will also prove useful to form 
charge and spin fields cj>i,p/<T — (0i,T ^ /V^- In terms 
of these variables. 



tf(lead) + # (dot 

2 ,0 



(3) 



^J2[dx (d x e lM f + (d x <j> i>a ) 2 

V / dx Ki (d x 9 hP ) 2 + — (d m <i 
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The Luttinger parameter in the reservoir, K%, is deter- 
mined by the interactions — K\ < 1 for repulsive inter- 
actions and K\ > 1 for attractive interactions. For the 
dot, K2 = 1. In this work, we will focus on the case of 
repulsive interactions, K± < 1. To analyze the physics it 
will prove useful to unfold the system, and work solely 
in terms of right moving fields. H Moreover, by forming 
linear combinations of the Bose fields in the dot and the 
reservoir, the system can be treated as two identical Lut- 
tinger liquids with an effective Luttinger parameter § 

K=™^. (4) 

The effects of Eq. | can be deduced by a renormaliza- 
tion group (RG) analysis. More generally, we will con- 
sider 

H int = t(r+ar^l s (0)^ s (0) + h.c) 

+ t'T^l^iM^M-hr*. (5) 




FIG. 2: RG flows of Eq. 



Though the t' term is not present in Eq. it will be 
generated upon renormalization. To lowest non-trivial 
order, the RG equations for the parameters are 



dXXV 1 Mfu 1 ^ 

r ]\h 

— = \ h - K(\ x y) 2 x h . 

dl 



(6) 



where X xy - t, X z - t', and X h - h. The RG flows in the 
X z — A^^-plane are plotted in Fig. ^. Notice that there is a 
critical point occurring when X xy = (1 — K)/{2yK). For 
yxy <- y^xy fa e coupling flows to zero, while for X xy > X xy 
the system flows to strong coupling. Since this critical 
point arises in the same way as in the pseudogap Kondo 
model, we will refer to it as the pseudogap Kondo critical 
point. M For X xy < X xy , the system flows to the fixed 
point where the dot is decoupled from the reservoir. (We 
will refer to this as the decoupled fixed point.) In terms of 
the effective Kondo model, the "impurity" is unscreened 
at low energies. For X xy > X xy , X xy initially decreases 
under the RG. However, it will eventually start to in- 
crease and then flow off to strong coupling. Integrating 
Eq. |, we find that X xy = 0(1) at a scale 



T K = E c exp 



-1 

w 



arccos(|<5|) — arctan(x /|^|)] 



(7) 



where x = (K - 1)/(2K) and \S\ = (X xy ) 2 /K — x\. 
For energies below Tk, the dot and the reservoir are 
strongly coupled. The strong coupling fixed point which 
arises is non-trivial — it corresponds to the 2-channel 
Kondo fixed point with a spin-1/2 impurity. Jll| This oc- 
curs because both spin-up and spin-down electrons in the 
reservoir try to occupy the single available charge state 
on the dot. 

It should be noted that a related system was considered 
recently in Ref. [IJ. In that work, the authors considered 
a resonant level coupled to a Luttinger liquid of spinless 
fermions. If the Luttinger parameter was smaller than 
some critical value, K < K C1 they too found a transition 
as one tuned the coupling between the dot and the Lut- 
tinger liquid. The authors of Ref. [l2| focussed on the zero 
temperature properties of their system. In this work, we 
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show that much rich physics can be observed at finite 
temperatures and frequencies. 

The quantity of experimental interest is the differen- 
tial capacitance. In terms of the effective Kondo model, 
this corresponds to the impurity susceptibility. || 
Hence, we will need to calculate correlation functions 
of impurity operators. To begin with, we will focus 
on the regime where X xy < 0(1). In this regime, 
we can calculate the impurity susceptibility using the 
RG. In general, an N-point impurity correlation func- 
tion Gjv(ti, • • • rjv; Aj, E c ) = (t z (ti) ■ ■ ■ t z {t n )) satisfies 
the RG equation 



G N (n,---T N ;Xi,E c ) = 0, (8) 



where Pi = dXi/dl, and 7 is the anomalous exponent. 
The solution of Eq. is 



G n (t u ■ ■ ■ t N ] A,, E c ) = exp 



-N I dl 7(0 




xG N (n,---T N ;Xi(l*),e 1 E c ). 
Using Eq. ^, we obtain 

G N (n, ■ ■ ■ r N] Aj, e c ) = e -"m-K)/2 e -NK*f ( n 



= 1*1 



where 

fin 
nn 

f(l*) = |*| tan 



xG N (T 1 ,---T N ;X i (l*),e- 1 E c ) , (9) 
+ {\S\ -x Q ) 2 e 2 ^ 1 ' 



y 2 -(\5\-x Q ) 2 e 2 \ s \ l * 
x 



1 - x I* 



\s\ r 



(\ xy = X xy ) 



arctan I 



1*1 



(X xy < X xy ) (10a) 
(10b) 

(X xy > X xy ) . 

(10c) 

In Eq. |j| x = iK- 1)/(2K ); y = X xy /VK and \5\ = 
y/xl -yi in Eq. [lOa} \S\ = y/(X xy ) 2 /K - xl in Eq. |0c . 
Choosing e l ~ E c /E, the correlation function on the 
right-hand-side of Eq. can be evaluated perturbatively. 

We start by considering the temperature dependence 
of the differential capacitance on resonance, C{h—* 0; T). 
Using Eq. || with N = 2, we obtain 

e K{K-l) 



C(h^0;T) 



-2K z f(T) 



AT 



(11) 



where f(T) is given by Eq. [u] with I* = \n(E c /ciT). (ci 
is an 0(1) constant.) From Eq. |ll|, we see that C(h — ► 
0;T) ~ T- 1 near the decoupled fixed point (T -> for 
X xy < X xy ). Moreover, near the pseudogap Kondo criti- 
cal point (X xy = X xy for T > E c exp[2K/(K-l)]) C(h -> 
0; T) ~ T At_1 where A T = (1 - K) 2 /2. It is also inter- 
esting to consider the differential capacitance at T = as 
a function of gate voltage, C(h;T = 0) = E c d(r z )/dh. 
Using Eq. ^ with N = 1, we obtain 



= K(Jf-l)/2 



±- 



-if 2 /(|ft|) 



(12) 
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FIG. 3: C{h;T = 0) vs. h for A a:!/ < X x c y (with if = 0.6). 
Solid line: \8\ = 0; dashed line: \5\ = 0.1; dash-dotted line: 
|5| = 0.2; dotted line: \S\ = 0.3. 



where the plus (minus) sign is for h > (/i < 0), and 
f(\h\) is given by Eq. with I* = \n{E c / c 2 \h\). (c 2 is 
another 0(1) constant.) The differential capacitance vs. 
gate voltage is plotted in Fig. [|. Near the decoupled fixed 
point, we find C{h\T = 0) ~ \h\ 2 \ s \~ 1 as |A| -> (|<5| is 
as in Eq. |l0a| ). Also, near the pseudogap Kondo critical 
point C(h;T = 0) ~ l/i^"" 1 where A h = (1 - if) 2 /4 
(H» £ c exp[2if/(Jf- 1)]). 

Notice that near the decoupled fixed point, C(h — > 
0; T) ~ T" 1 as T -> 0. This Curie- Weiss-like form arises 
because the "impurity" behaves basically like a free spin. 
However, the local "moment" is reduced from its non- 
interacting value. 0, [l2| [l3| From Eqs. [l2[ we see that 
the amount by which the local "moment" is reduced de- 
pends on K, as well as the value of X xy . Also, notice that 
Ay = 2A/j near the pseudogap Kondo critical point. In 
Ref. it was shown that the exponents near the criti- 
cal point of the pseudogap Kondo model satisfy certain 
hyperscaling relations. To the order of accuracy that we 
have worked, our results are consistent with these hyper- 
scaling relations. 

Now, we consider the physics in the regime X xy > X xy 
for energies below Tk (with Tk given by Eq. fj]). In 
this regime, the system is close to the 2-channel Kondo 
fixed point. To proceed, we follow Ref. |l4| and form 
combinations of the fields in the dot and the reservoir: 
4>R,c 4>R,s P , 4>R,f, and (j>R tS f- Then, we perform the 
unitary transformation, U = exp (i^Air/K T z (f)R t fi0) \ , 
which ties charge to the "impurity" . Finally, we intro- 
duce new fermion fields, d ~ r _ , X ~ e ^V^<t>R, S f ^ anc j 
/ ~ e l ^f^4>n,s _ Upon performing these transformations, 
i?i n t becomes 

flint = v F X xy id^+d) (A+(0) - X(0)) 

+ v f ^/k(x z -l) (Sd -l/2)/t(0)/(0) 

- v F X h idU - 1/2) , (13) 

where X xy , X z , and X h are the renormalized values of the 
couplings. 
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O < < • ) ) • < 

decoupled pseudogap 2-channel 
Kondo Kondo 

FIG. 4: Fixed points and RG flows for the pseudogap Kondo 
model. 

Using Eq. |l3|, we can calculate the differential capac- 
itance near the 2-channel Kondo fixed point. Starting 
with the differential capacitance on resonance, we find 
(ignoring the irrelevant (A 2 — 1) term) 

C(h^0;u,T)= (14) 

J_ f ^tanh ( xTk \ _J 1 

T K J 2ir V 2T ) x 2 + 1 x - (u/T K ) - i0+ ' 

For uj — and T <C Tk, this reduces to the well-known 
result for the impurity susceptibility of the 2-channel 
Kondo model C(h -> 0;T) = l/(7rT Kr ) ln(T K /T). We 
can also calculate C(h;T — 0). Using Eq. |l3| (ignoring 
the irrelevant (A z — 1) term) 



h 




r xT K \ 




^ — tanh | 




2vr * 


v 2T J 



(i 2 -(V^) 2 ) 2 + i 2 ' 



(15) 

For T ~ and \h\ < T K , C(/i;T = 0) = 
4:/(nT K )ln(T K /\h\). Notice that C(/i -> 0;T) (C(/i;T = 
0)) diverges as T — > — * 0). However, the divergence 
in this case is weaker than what occurs near the decou- 
pled fixed point. This is because, near the 2-channel 
Kondo fixed point, charge is tied to the "impurity". As 
a result, the ground states r z = 1/2 and r z = —1/2 
are orthogonal, in that they are not connected by t + or 
r~ . [ fEB] This removes the power- law divergence which oc- 
curs near the decoupled fixed point, and replaces it with 
the weaker logarithmic divergence. 

In the above discussion, we saw three fixed points arise 
(shown schematically in Fig. ||): (1) the decoupled fixed 
point, (2) the pseudogap Kondo critical point, and (3) 
the 2-channel Kondo fixed point. The pseudogap Kondo 
critical point and the 2-channel Kondo fixed points are 



particularly interesting because they are non-trivial scale 
invariant fixed points. As a consequence, one should be 
able to observe w/T-scaling near these fixed points by 
applying an AC component to the gate voltage. More 
specifically, we expect the dynamical capacitance on res- 
onance to have the form 

C{h^Q]U),T) = T v - 1 X{u)/T). (16) 

Near the pseudogap Kondo critical point, we can calcu- 
late the scaling function X(ui/T) for (1 — K) < 1, In 
the leading logarithm approximation, we find v = At 
(A t = (1 - K) 2 /2) and 



X(u>/T) = 




where B is the beta function. Near the 2-channel Kondo 
fixed point, we use Eq. |l4] to obtain v — 1 and 

X ^ = ln (n^T)) + 2k tanh (£) ' 

Note that X(co/T) is, in general, complex. Therefore, 
the differential capacitance will have components both 
in-phase and out-of-phase with the gate voltage. 

To summarize, a (large) quantum dot coupled to an 
interacting one-dimensional electron liquid could provide 
a controlled environment to observe pseudogap Kondo 
physics. By tuning the coupling between the dot and the 
one-dimensional electron liquid, one can access the vari- 
ous fixed points which arise: the decoupled fixed point, 
the pseudogap Kondo critical point, and the 2-channcl 
Kondo fixed point. Moreover, this system provides the 
remarkable opportunity to directly probe impurity prop- 
erties via differential capacitance measurements. As the 
differential capacitance of a large quantum dot has re- 
cently been measured, [|l6) we are hopeful that the physics 
described in this work can be observed in the near future. 
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